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We  Jer'vc  ihe  umc-deix*ndcM  spectra  of  a  laser-driven,  imeravtmg.  collection  of  spauall>  distributed  two-level  atoms  We 
assume  a  s>siem  with  dimensions  large  relative  to  a  wavelength  uhieh  is  imtiallv  in  the  alomie  and  electromagnetic  field  ground 
state  and  at  the  initial  time  begins  to  interact  with  the  coherent,  cxtemalb  applied  radiation  field.  The  explicit  spaeial  disti  thuiton 
of  the  ivs ts-lev el  atoms  affects  the  lime-dependent  spectra  of  the  transmitted  radiation  and  the  manybodv  interaction  leads  to  the 
appearance  of  additional  sidebands 


The  fluorescence  spectrum  of  a  system  of  a  large 
number  of  tv.o-lev el  atoms  driven  b>  a  strong  reso¬ 
nant  field  ha.  been  treated  for  atoms  that  are  con¬ 
tained  in  a  region  that  is  small  relative  to  a 
wavelength.  For  high  laser  intensities,  resonances 
displaced  from  the  line  center  by  twice  the  Rabi  fre¬ 
quency  were  predicted  ( 1-3)  in  addition  to  the  usual 
Stark  sidebands  [4.5],  The  intensity  of  the  addi¬ 
tional  sidebands  in  resonance  fluorescence  experi¬ 
ments  is  expected  to  be  very  weak  and.  so  far.  they 
have  not  been  observed. 

In  this  communication,  we  treat  the  appearance  of 
additional  sidebands  in  a  time-dependent  spectrum 
[6-9]  for  a  nonstationary  state  of  the  two-level  sys¬ 
tem.  We  follow  here  a  previous  model  [  10]  which  is 
based  upon  a  quantum  electrodynamical  represen¬ 
tation  of  a  spaeial  distribution  of  a  large  number  of 
two-level  atoms  interacting  with  the  quantized  elec¬ 
tromagnetic  field.  In  the  previous  work  [10]  we 
studied  the  steady-state  conditions  for  this  system. 
We  picture  that  different  individual  atoms  follow 
Rabi  oscillations  with  the  same  frequency  but  differ¬ 
ent  phases  so  that  by  ensemble-averaging  we  get  many 
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atom  steady  states.  In  the  present  work  we  assume  a 
sample  of  two-level  atoms  with  dimensions  very  large 
relative  to  a  wavelength  which  is  initially  in  the 
ground  state  and  at  time  /  -  0  begins  to  interact  with 
the  coherent  externally  applied  resonant  radiation 
field.  We  studv  the  spectrum  of  the  radiation  trans¬ 
mitted  through  the  sample  of  two-level  atoms  and 
show  that  the  explicit  spaeial  distribution  of  the 
atoms  affects  the  shape  of  the  radiation  spectrum  and 
can  lead  to  the  appearance  of  additional  sidebands 
[II].  The  additional  sidebands  should  follow  from 
modulation  of  the  collective  dipole  oscillations  at 
frequency.  2 1  |  in  addition  to  modulation  at  fre¬ 

quency  j  ioK  |  which  leads  to  the  fundamental  bands. 

We  simplify  the  present  treatment  for  time-depend¬ 
ent  spectra  by  assuming  that  the  filter  time  [  6—9  j  is 
short  relative  to  the  relaxation  and  dephasing  times, 
but  long  relative  to  the  period  time  of  Rabi  oscilla¬ 
tions.  finder  these  conditions  we  get  a  quasistation¬ 
ary  spectrum  which  is  damped  slowly  in  time.  In  the 
present  approach  we  describe  the  time-dependent 
collective  dipole  oscillations  by  using  a  semiclassical 
model  and  treat  the  nonlinear  terms  by  a  perturba¬ 
tive  approach  in  which  the  reaction  field  is  taken  as 
small  relative  to  the  externally  applied  field.  Follow¬ 
ing  the  above  considerations,  the  present  analysis  is 
good  as  long  as  the  magnitude  of  the  collective  oscil- 
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laiions  is  large  relative  to  other  sources  of  fluctua¬ 
tions  including  those  related  to  the  Langevin  terms 
that  represent  quantum  fluctuations. 

The  hamiltonian  which  describes  the  system  is 
given  in  the  rotating  wave  and  in  the  electric  dipole- 
approximation  by 

//=//„  +  //  . 

\ 

W,,  =  4ftw  y  a'/' +fi  V  OJ1.C1I Ui  . 

1  I  A 

*/' =  -t£  I  exp(  iAr  ) 

a.  A  1  ! 

\h  K 

-  —  X  M k ff' • '  exP[  -( -*-r, )]  +h.c. .  ( I ) 

where  h.c.  denotes  hermitian  conjugate  and 
g[‘'  =(2nojk/fil’)'  -Pp,-i 

//,,  includes  the  free  atoms  and  free  Held  hamiltoni- 
ans.  We  use  here  the  usual  operators  related  to  the 
SU(2)  algebra  fora  single  atom:  <r"‘  represents  the 
inversion  of  population  operator  and  o':'  are  the 
raising  and  lowering  operators  for  atom  i  with  coor¬ 
dinates  r,.  is  the  Rabi  rate  associated  with  the 
applied  coherent  Held  amplitude,  to.,  and  k  are  the 
carrier  frequency  and  wavevector  of  the  applied  field, 
respectively.  I  is  the  quantization  volume  for  the 
field,  w,  is  the  frequency  of  the  mode  A.  and  x  is  a 
unit  vector  in  the  direction  of  polarization  of  the 
applied  radiation  field  that  is  being  taken  to  be  in  a 
coherent  state  and  propagating  along  the  r-axis.  /’  is 
the  dipole  matrix  element  and  p,  is  the  unit  vector  for 
the  dipole  of  atom  /. 

In  the  previous  work  [10]  we  developed  the  Hei¬ 
senberg  equations  of  motion  in  the  bad  cavity  limit 
for  this  system.  The  rapid  time  and  spacial  depen¬ 
dence  of  the  dipole  operators  has  been  eliminated  by 
the  transformation: 

ol:\,(l)=o'l'(i)  exp[  — i (w„/-*„-r,)]  .  (2) 

By  using  the  coupled  hierarchy  of  the  Heisenberg 
equations  of  motions,  and  by  adiabatically  eliminat¬ 
ing  the  variables  associated  with  the  field  modes  one 
gets  terms  in  the  equations  of  “spontaneous",  “coop¬ 
erative".  “induced",  and  “Langevin  force"  origin. 
While  the  “spontaneous”,  "induced"  and  "Langevin 


force"  terms  depend  only  on  single  atom  response, 
the  "cooperative"  terms  include  factors  of  operators 
belonging  to  pairs  of  different  atoms  with  arguments 
expressive  of  retardation.  The  basic  equations  are 
obtained  by  taking  the  expectation  values  and  facto¬ 
rizing  the  multiplication  of  the  dipole  operators. 
Factorization  is  a  justifiable  approximation  in  a 
many-body  theory  as  the  correlation  between  each 
pair  of  atoms  tends  to  be  overwhelmed  by  the  total 
interaction  of  each  atom  with  the  many-atom  sys¬ 
tem.  The  basic  equations  for  a  long  sample  of  two- 
level  atoms  with  dimensions  larger  than  a  wave¬ 
length  are  obtained  [10]  by  describing  the  expecta¬ 
tion  values  of  the  atomic  operators  as  continuous 
variables  and  defining  < cr  . r.  f )  >  and  < r Mr.  t)  > . 
respectively,  as  the  complex  dipole  and  the  inversion 
of  population  per  unit  volume  at  point  r.  Follow  ing 
the  previous  work  [10]  we  get  the  following  equations 

dv  o  (r.  /)  >/d/=  -/y ,  [  <ff  (r.  t)  >  +/i] 

- ;  ( 20*(  r.n  +  cju  ]  <  <r  .  „(  r.  t ) )  +  c.c !  .  (31 

d(<; .  „(r.  /)  >/dr=tiJ-/(;)<ff  .  „(r.  I)  > 

+  [f3(r.  O+wR/2] <<7. (»■./)>  ,  (4) 

where  0(r.  i)  is  a  continuous  variable  corresponding 
to  the  reaction  field  and  is  given  by  [  10] 

H(r.  t)  =(  3/t/y/A  'i 

x  |  dr'  <cr . „(r  . (r-r  )/<•)  >  .  (5) 

where  n  is  the  number  of  two-level  atoms  per  unit 
volume  and  .)  =  «;  -o,,  is  the  deviation  of  the  applied 
field  frequency  from  resonance.  Self-consistently 
from  the  quantum  model  /(,  =/i.  /i  .-/II2  where  / i  is 
the  spontaneous  decay  constant.  We  may  introduce 
additional  homogeneous  broadening  by  considering 
IU  and  /F  as  empirical  constants.  Eqs.  (3)^(5)  lead 
tc  a  set  of  equations  which  arc  of  the  form  of  the  usual 
Maxwell  Bloch  equations  if  we  perform  the  space 
derivative  in  eq.  i  5 )  and  neglect  the  renormalization 
of  the  resonance  frequency  caused  by  coherent  dipole- 
dipole  interaction  [10]. 

We  calculate  the  time  dependence  of  the  total  elec¬ 
tromagnetic  field  transmitted  at  the  end  of  a  short 
sample  for  which  the  conditions 
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1 0(r,  r)l  <$  |wKl.  Pi  :<  Iwk  )  • 

are  fulfilled.  Under  these  conditions  the  reaction  field 
is  given  approximately  by 

0(r,  I )  =  «<(7 . t)  >  •  (6) 

where 

o;  =  (3  nfi/k'):.  (7) 

In  order  to  explain  further  the  use  of  eq.  (6).  we 
substitute  this  equation  in  the  ordinary  Maxwell 
equation: 

c  '  r)ti( i)/< + 

=  ( t)  >  (8) 

and  find  that  eq.  ( 6 )  is  valid  it 

( f/ ( (7  ()>/('/;),  —  0  .  (9a) 

or  more  explicitly  if  consistent  with  our  first  order 
perturbation  theory: 

|  r[  it  logv <7  . >/ifr],  I  <  I  (9b) 

where  i  = t-:!c  is  the  retarded  time.  Since  vc  is  a 
small  parameter,  the  derivative  of  <<7  .  „>  relative  to 
-  will  contribute  only  to  second  order  corrections 
which  can  be  neglected.  The  restriction  of  the  discus¬ 
sion  to  a  strong  field  in  which  our  approximations 
are  valid  simplifies  the  present  treatment  for  the 
appearance  of  the  additional  sidebands. 

By  assuming  exact  resonance  ( _1  =  0)  and  substi¬ 
tuting  eq.  (6)  into  eqs.  (3)  and  (4).  we  get 

<cr  . )  =  —/?,[  (a  >  +n)  -  Wk<<7  .  „> 

-to *<<7  „> -4«<cr  .„>  <<t  d>  ,  (10) 

<(7.„>=-/U<(7.„>+itoj$<f7r> 


in  a  and  (i,  2-  We  get 

<*_->  +  |Wr  I :  <a:  >  - 3«<ff_-  >  0-.  > 

+  (P,  +/?,)<  <7- >=0.  02) 

with  the  first-order  solution 

<<7,(r)>  =<7-0  cos(  |toK  I?)  expj  —  [(ft,  +0:/2)]f  } 

+  y ( <v/ 1 to k  l)<7-„  sih(2|wR  |/)  exp[  -  ((i,  +/MC]  ■ 

(13) 

We  have  assumed  here  ’hat  <er.-(x)>=0,  and 
<<7  (  x  )>  =0  for  a  strong  fidd  at  exact  resonance. 

Eqs.  (12)  and  (13)  derived  here  are  similar  to  the 
equations  derived  by  Senitzky  [  1  ].  The  assumption 
of  conservation  of  total  angular  momentum  made  by 
him  is  not  valid  in  our  treatment  due  to  the  damping 
mechanism.  The  parameter  a  derived  here  is  for  a 
sample  with  dimensions  large  relative  to  a  wave¬ 
length.  It  depends  on  r  via  eq.  ( 7 )  and  is  essentially 
different  from  a  corresponding  parameter  a  derived 
in  ref.  [  1  ]  for  a  sample  with  dimensions  smaller  than 
a  wavelength.  The  broadening  mechanism  which  was 
neglected  there  [  1  ]  is  introduced  here  via  the  param¬ 
eter  p,  +fi:. 

The  solution  for  {a  .„(/)>  of  first  order  is  given 
by 

<<7.,. >  =  I w*  a,„  expj  -  HP,  +/}:)/2]f) 

sin(  |  toR  |/')  jg  i  1  -cos(  |  toK  If  A 
I  ton  I  2  |wKi;  1 

+  or<j:„(toj*/2)  exp[  -/?,  +/U)/'] 

1 1  —  cos(  2 1  toR  |f) 

\  2  |wK  | : 


+  «<(7  .  (11) 

where  the  dot  represents  the  time  derivative  d/dt 
Here  we  have  used  a  short  notation  in  which  it  is 
undeistuod  that  all  expectation  values  are  functions 
of  retarded  time  I'  and  at  position  r.  The  explicit 
dependence  of  the  expectation  values  on  r  enters  only 
via  the  parameter  a. 

Assuming  /?,  >«  |toK|.  a««|wK|.  we  solve  eqs. 
(10)  and  (I  I )  by  iterative  substitutions  in  which  we 
neglect  all  terms  which  are  of  scconJ  order  or  higher 


According  to  eqs.  (13)  and  (14)  the  oscillation  is 
modulated  at  frequency  2|wR|  in  addition  to  the 
modulation  at  frequency  |wK|. 

Assuming  that  the  CW  laser  is  activated  at  t'^0 
and  that  the  two-level  atoms  are  initially  in  the  lowrr 
state,  the  time-dependent  spectrum  is  given  by  [7] 
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S=4y  Re  j  dr  e\p[(;’-Llco)r] 

O 

x  j  do  exp[ -2;'(/-o)] 


X  <  <£*(/  +T)  £(/;)>>  . 

(15) 

where,  for  our  case 

£(f )  =  [««  -t-20*(;.  I)]fi/P . 

(16) 

and  My  is  the  filter  time. 

By  assuming  the  conditions 


/?,  .  (17) 

the  ensemble  average  in  eq.  (15)  refers  to  averaging 
over  the  Rabi  oscillations  where  the  electromagnetic 
field  at  time /'  is  obtained  by  using  eqs.  ( 6 ).  (14). 
and  (16).  We  notice  that  the  zeroth  order  and  first- 
order  value  for  the  complex  dipole  introduce,  respec¬ 
tively.  first-order  and  second-order  corrections  to  the 
zeroth  order  electromagnetic  field  which  is  propor¬ 
tional  to  wK.  We  can  distinguish  in  the  expression  for 
the  electromagnetic  field  between  terms  that  are  con¬ 
stants.  those  that  oscillate  ai  frequency  |w„|  and 
those  that  oscillate  at  frequency  2|mr|.  Due  to  the 
ensemble  averaging  over  the  Rabi  oscillations  the 
terms  in  the  correlation  function  which  couple  dif¬ 
ferent  oscillation  frequencies  can  be  neglected. 

In  order  to  simplify  the  expression  for  the  correla¬ 
tion  function  we  include  in  each  term  that  oscillates 
at  a  given  frequency  only  that  part  which  results  from 
the  lowest  order  in  perturbation  theory.  Therefore  we 
neglect  the  first-  and  second-order  corrections  to  the 
constant  term  wR  and  we  neglect  the  second-order 
correction  to  the  terms  modulated  at  frequency  |oiR  | . 
Following  these  considerations  we  obtain 

<<£(/:  +T)*E(t2))}  =  (>r/P:)\wR  |- 

X  { 1  +ji<5:  cos(|wR|r)  exp[  -  (0,  +p2)(t:  +  r/2)] 

+  exp[  -(/?,  +P:){2l2  +  x)]  cos(2|wr  |t)}  , 

(18) 

where 

<5  =  «<».„/ 1 aiR  |  .  (19) 


In  evaluating  the  spectrum  according  to  eqs.  (15) 
and  (18)  we  use  the  approximation  /?,;«'/  and 
neglect  terms  which  include  exponents  with  argu¬ 
ments  as  products  of  r  and  /?,  or  but  take  into 
account  terms  which  include  exponents  having  argu¬ 
ments  as  multiplications  of  t  and  /f,  or  as  they 
will  affect  the  decay  of  the  spectra  for  long  limes  t  . 
Also  we  develop  the  expression  for  times  t  >  !/;■  so 
that  exp(  -  2)7 ' )  ->0.  By  following  these  approxima¬ 
tions  [7]  we  obtain  by  a  straightforward  calculation 

c;  2|wr  {'fi1  \  y 

P:  l(a>  -to,,)2  +;,: 

+  l4S:  exp[  -  (/? ,  +P:)l  } 


(<o-a>,,+  |wR  | ) -  +;,: 


Uo-mu-  |wR  |  )•’+;•'  J 

+  exp[  -2 (fi,  +p:)f] 


(  cu  -  co„  +  2|ozr  | ):  +  y ' 


-r  -  .  .  if  -  \  ; 

Uo  -  CD,,  -2](Dr  |  )■  +y  J  ) 

Eq.  (20)  shows  that  within  our  approximations  the 
time-dependent  spectra  for  the  transmitted  radia¬ 
tion  includes  the  laser  line,  and  the  fundamental  and 
the  additional  sidebands  with  lorentzian  line  shapes 
whose  centers  are  located,  respectively,  at 
co,i±  |<wK|  and  <y„±2|wR|.  Since  according  to  our 
conditions  y  dominates  the  correlation  function,  the 
width  of  each  band  is  fixed  by  the  filter  time  l/;1  and 
is  different  from  the  Mollow  line  shape  in  which  the 
width  of  each  line  is  related  to  /?.  The  lorentzian  fun¬ 
damental  and  additional  sidebands  are  decaying, 
respectively,  with  the  rates  pt  +f}2  and  2(0 ,+P:).  By 
our  perturbative  theory  we  find  that  the  ratio  of  the 
intensity  of  the  fundamental  sideband  to  the  mono¬ 
chromatic  laser  line,  and  that  of  the  additional  side¬ 
band  relative  to  the  fundamental  is  given  by  S2/ 4 
(similar  to  a  relation  found  in  ref.  [11). 
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one  gets  shifted  energy  levels  with  the  corresponding 
frequencies  ojn.  w„±|<uR|.  The  additional  side¬ 
bands  result  from  transitions  in  which  the  dipole 
moment  is  doubled.  In  order  that  such  transitions 
occur,  the  dipoles  of  different  atoms  must  have  a  non- 
negligible  probability  to  oscillate  in  phase  with  one 
another.  While  the  probability  for  this  kind  of  tran¬ 
sition  is  very  small  for  the  many -atom  steady  states, 
it  becomes  relatively  strong  in  the  case  treated  here 
due  to  the  transient  cooperation  induced  by  our  ini¬ 
tial  conditions.  For  long  times  the  system  tends  to  a 
steady  state  for  which  the  relative  phases  of  different 
dipoles  become  randomized.  Therefore  we  suggest 
that  the  observation  of  additional  sidebands  can  be 
made  in  a  time-dependent  spectrum. 
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